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There are different techniques to sense the wavefront phase-distortions due 
to atmospheric turbulence. Curvature sensors are practical in their sensitivity 
being adjustable to the prevailing atmospheric conditions. Even at the best 
sites, the turbulence intensity has been found to vary at times over only a 
few minutes and regularly over longer periods. Two methods to automatically 
adjust the sensitivity of a curvature sensor are proposed: First, the defocus 
distance can be adjusted prior to the adaptive-optics (AO) loop through the 
acquisition of a long exposure image and can then be kept constant. Secondly, 
the defocus distance can be changed during the AO loop, based on the voltage 
values sent to the deformable mirror. We demonstrate that the performance 
increase - assessed in terms of the image Strehl-ratio - can be significant. @ 
2010 Optical Society of America 
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1. Introduction 

Curvature wavefront sensors [1] have a number of advantages over the more commonly 
used Shack-Hartmann sensors [2]. One useful characteristic that has not often been used in 
practice is their adjustable sensitivity to phase distortions: While Shack-Hartmann sensors 
are designed for one fluctuation amplitude taken to be prevalent during the instrument's 
lifetime, the dynamic range of curvature sensors varies with the distance of the detector 
from the focal plane (the defocus distance). Accordingly, curvature sensors can be tuned to be 
optimally sensitive to the wavefront aberrations that prevail during a particular observation. 

In most adaptive-optics applications the amplitude of the phase fluctuations does indeed 
vary substantially. In astronomy, the seeing - a measure of the phase-fluctuation amplitudes 
- varies on a time scale of typically several minutes, which is shorter than the integration 
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time of most astronomical observations. The seeing varies even faster and might become 
prohibitive when an observed object, such as a satellite, moves across the sky faster than 
astronomical targets. 

In 1991, F. Roddier [3] suggested to use a feedback loop to keep the rms tip-tilt signal-error 
constant. The sensor would then automatically adjust itself to changing seeing conditions 
and could always operate in its optimum linear range. Existing curvature-sensor based AO 
systems do not, as yet, implement such a feedback loop. For example on PUEO, the 19 
element curvature sensor installed behind the Canda-France-Hawaii telescope, the defocus 
distance (termed optical gain in the PUEO manual) is not, generally, adjusted to changing 
atmospheric conditions: "The optical gain value is (...) by default set to 128 for a point-like 
object: this value may need to be reduced for an extended object or a double star" [4]. 
On MACAO [5], the curvature AO systems installed behind the Very Large Telescope, the 
defocus distance is adjusted as the reference source is changed. The adjustment depends on 
the seeing measured by the DIMM-monitor [6] and on the apparent size of the reference 
source (personal communication). 

We propose two practical implementations of defocus feedback-loops: First the defocus 
distance can be adjusted prior to the AO loop through the acquisition of a long exposure 
image and can then be kept constant. This solution is similar to the MACAO approach. 
Secondly, the defocus distance can be changed during the AO loop, based on the voltage 
values sent to the deformable mirror. We analyze the resultant performance increase through 
numerical simulations. 

2. Optimal defocus distance 

The output signal of a curvature wavefront sensor is the fractional intensity difference v = 
(x — y)j (x + y) , where x and y are the photon numbers registered on either side of the focal 
plane. In closed-loop AO observations v tends to be near and can then be expressed as a 
function of the wavefront phase, [7]: 

v(f) = " A ^" [V(P(fl) • V(<Kfl) + P( jO " V 2 0(f-)] (1) 

f is the position vector in the detector plane. / is the focal length, I the distance between 
the detector and the focal plane, and f /If is thus the position vector in the pupil plane. 
As is the sensing wavelength, and P - the pupil transmission function - equals 1 inside and 
outside the pupil. Eq. 1 makes use of the Fresnel approximation (i.e. D 2 /((f — l)X) > 1) 
and it assumes the size of the diffraction patterns to be negligible compared to the magni- 
tude of the intensity fluctuations in the detector plane. This condition, the geometric-optics 
approximation, constrains the defocus distance to As/ 2 « r o(^s)l- The Fried parame- 
ter, ro(Ag), characterizes the turbulence intensity in dependence on the wavelength, As. In 
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regions where the pupil transmission is uniform, the measured signal is the Laplacian of 
the wavefront phase. On the pupil edge - i.a. along the central obstruction and along the 
telescope spiders - curvature sensors measure curvature and radial tilt. 

In line with Eq. 1 we want small defocus distances, I, in order to have sufficient sensitiv- 
ity to small phase fluctuations, but for a given intensity spatial fluctuation size within the 
aperture, d = D/yN , we want / sufficiently large in order to have dl/f larger than the 
diffraction pattern \s(f — 0/ r o(^s)- As Rigaut et al. [8] state, to limit the aliasing effect 
the diffraction should be large to filter out the not measured high-order modes of the tur- 
bulence. Accordingly the optimal distance is dl/f — k\ \s{f — 0/ r o(^s)- The coefficient k\ 
is sufficiently small to minimize aliasing effects and sufficiently large to avoid non linearities 
due to diffraction. Let Iq be the optimal defocus distance. The defocus length is typically 
substantially smaller than the telescope focal length, I << /, so that: 



the coefficient k\ depends on the exact size of the sub-apertures, and in particular on the 
size of the central obstruction. But, for a given AO system, k\ is independent of the observing 
conditions such as seeing and stellar magnitude. 

3. Defocus adjustment prior to the AO loop 

Here and in the subsequent section, numerical simulations of curvature AO systems will be 
presented. In this section the seeing is kept constant during each sequence and the point is 
explored, whether the optimal defocus distance can be determined from the diameter of a 
long exposure image recorded prior to the AO loop. 

The simulated curvature AO-system has N = 80 sub-apertures, which is an aver- 
age of existing astronomical curvature systems (PUEO/Canada- France-Hawaii telescope: 
19 [9], MACAO /Very Large Telescope: 60 [5], NICI/Gemini-south telescope: 85 [10], AO- 
188/Subaru telescope: 188 [11]). The sub-apertures are placed on concentric rings as shown 
on Fig. 1, and cover equal fractions of the pupil-surface. 

The numerical simulations are based on the code developed by F. Rigaut 
(www.maumae.net/yao/). For 22 seeing values between 0.2" and 1.4" (at A = 0.5 /im) a 
matrix of incoming phase values (the phase screen) is generated and is shifted stepwise in 
front of the pupil. The distribution of phase values follows the Kolmogorov model. The 
step-size per iteration is indicated in Table 1. Observation sequences are then simulated over 
3 seconds. 

Initially, the optimal defocus length, I, and the loop gain, g, are determined by repeating 
the numerical simulations for g = 0.1,0.2, ...1 and I = 0.1,0.2, ...1.8 m in an f/QO beam. A 
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separate command matrix, which converts a set of N sensor measurements into N voltage 
values, is computed for each defocus distance. These matrices are computed by first generat- 
ing N influence functions: N mirror surfaces corresponding to a unit voltage applied to each 
of the actuators. The influence functions are independent of the defocus distance and are 
generated only once. Note that the influence functions are stored in a 237 x 237 matrix, while 
the pupil-diameter equals 220 pixels: all points of the influence functions that lie outside the 
pupil are set to zero. For each defocus distance, Yao computes an interaction matrix through 
Fresnel diffraction: the N x N matrix contains the N sensor measurements induced by a unit 
voltage applied to each of the N actuators. The command matrix is obtained by inverting 
the interaction matrix. 

Fig. 2 shows the Strehl ratios attained at 1.65 /xm in dependence on the seeing value at 
0.5 /xm and compares them to the ratios achieved when I is kept constant while the seeing 
changes. In the latter case the loop gain is still optimized for each seeing value. As shown 
in Fig. 2, the correction performance will be especially low if the curvature sensor is tuned 
for excellent atmospheric conditions (i.e. the defocus distance is small). The reverse effect - 
a performance drop due to an overly pessimistic sensor adjustment - also exists, but is less 
fatal: for excellent seeing conditions, 8q = 0.2", the Strehl ratio decreases from 0.94 (7=0.2 m) 
to 0.87 (7=0.6 m) and 0.83 (7=1.1 m). 

Since it is impractical to determine the best defocus distance, Iq, by running a series of 
AO loops during actual observations, we need to explore whether - from the diameter of a 
long-exposure image acquired prior to the AO loop - Z can be directly obtained through the 
formula: 



As and A/ are the sensing and imaging wavelengths. 9 50 (\j) ~ A//r (A/) [rad] is the diameter 
of the circle that contains 50% of the intensity on a long-exposure image at wavelength A/. 
The coefficient k 2 depends on instrumental parameters such as the central pupil obstruction 
and it needs to be determined through simulations or instrumental calibrations. D is the 
pupil diameter, / the telescope focal length, and N the number of sub-apertures. Table 1 
gives their numerical values. 

The coefficient k 2 is determined via Eq. 3: 22 values of Iq and 6*50 are obtained for seeing 
values between 0.2" and 1.4", at A = 0.5 /xm. A linear fit yields: k 2 = 0.68 ± 0.07. As 
seen on Fig. 2, the performance of the AO loop is then the same - whether the defocus is 
determined prior to the AO sequence from the diameter of a long-exposure image (solid 
line), or whether many different defocus values are tested (triangles). It is thus possible 
to optimize the system's performance by linking the defocus distance to the diameter of 
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a previously recorded, long exposure image. Another possibility is to adjust the defocus 
distance in dependence on the seeing measured by a turbulence monitor. On MACAO [5] 
the defocus is adjusted in terms of the seeing as measured by the DIMM [6]. The disadvantage 
of this approach is, that the seeing predicted by the DIMM monitor is often much worse than 
the seeing sensed by the Very Large Telescopes [12]. This seems to be explained by a highly 
turbulent layer close to the ground, which is seen by the DIMM (installed on a 6m high 
tower), while it does not affect the telescopes (primary mirrors at 11m above the ground, 
top of enclosures at 30 m). 

If the magnitude of the wavefront aberrations changes significantly over the course of the 
observation, the defocus length needs to be adjusted during the AO loop. In the following 
section a suitable procedure is suggested and is tested via numerical simulations. 

4. Adjustment during the AO loop 

Fig. 3 shows seeing values recorded during two nights by the DIMM [6] and MASS [13] 
instruments at the Mauna Kea observatory. Such significant and rapid seeing variations have 
also been reported at the Paranal observatory by Rigaut and Sarazin [14]. The significant 
and rapid seeing variations suggest that even observations with integration times of just 
several minutes would benefit from a defocus feedback loop. 

4. A. Method 

In the previous section, the defocus distance has been assumed to be determined prior to the 
AO loop from the diameter of a long-exposure image. To adjust it during the observation, 
a constantly updated estimate of the seeing is required. This estimate is here taken to be 
obtained from the voltage values sent to the deformable mirror: 

The phase of the wave-front is continuously redetermined from the voltages and the influ- 
ence functions of the actuators: = 21 V. The voltage vector V has N elements, N being 
the number of actuators. The matrix I [rad/V] contains the N influence functions. Each 
influence function is sampled over M x M points, and I is therefore of size: M x M x N. <fi 
[rad] is the wavefront phase at wavelength A5 and is likewise sampled over M x M points. 

The phase variance, averaged over a certain number of loop cycles, K, is then used to infer 
the Fried parameter: 

1 1 K M ( D \ 5/3 

k=l m,n=l v uv OJ 7 

a\ is the average phase variance over K cycles at wavelength Xg. <pk is the wavefront phase 
at wavelength A5 for cycle number k. The exact value of k need not be assessed, since the 
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defocus distance is related to D/r (Xs) through another coefficient (see Eq. 2). So that: 



&3 is then determined from simulations or instrumental calibrations. The actual defocus 
distance is set to the closest value for which a command matrix has been computed. 

4-B. Numerical simulations 

At the outset, 30 command matrices are obtained for defocus values in the [0.1, 3.0] m range, 
with a 0.1m increment. The computation of the command matrices follows the principle 
described in Sec. 3. An initial defocus value is specified in terms of Eq. 3 from a long exposure 
image acquired with the initial seeing value. 

The coefficient k 3 is obtained via Eq. 5: Twenty-two AO sequences are simulated for seeing 
values between 0.2" and 1.4", at A = 0.5 /im. Each sequence contains 500 iterations. The 
optimal defocus distances corresponding to those seeing values have already been determined 
in Section 3. Starting at the 300 th iteration, the average phase variance over the last 100 
iterations - i.e. during the last 0.5 s - is computed. The mean and standard deviation of 
the 200 realizations is determined, and the coefficient, k 3 , is obtained through a linear fit: 
k 3 = 1.2 ±0.1. 

During the AO loop the phase screen is shifted stepwise in front of the pupil. The size of 
the matrix and the step-size per iteration are indicated in Table 1, and allow for a maximum 
exposure-time of 70 seconds. Every 5 seconds the values of the phase screen are scaled to a 
new seeing value. In order to adjust the defocus distance during the AO loop, the wavefront 
phase variance is computed after each iteration from the eighty voltage-values and their 
associated influence functions. 

Figures 4 and 5 compare the performance of the AO system for two approaches: 

1. The defocus distance is kept constant. This is the approach presented in Section 3, but 
in the present treatment the seeing varies during the observing sequence. 

2. Starting from the 300 th iteration, the defocus distance is updated every fourth iteration: 
the average phase variance over the last 100 iterations - i.e. during the last 0.5 s - is 
used to determine a new defocus distance through Eq. 5. The defocus is set to the 
closest value for which an interaction matrix has been computed, i.e. it is set to 0.1, 
0.2,... 2.9 or 3.0 m. 

In both cases, the AO loop gain is adjusted for the average seeing over the sequence and 
is then kept constant. In Fig. 4 the seeing increases continuously from 0.2" to 1.4". The 
long exposure image is acquired when 9 = 0.2" and the initial defocus length is thus small 
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(0.2 m). Under this condition the Strehl value breaks down when the seeing increases while 
the defocus distance is kept constant at its low value. Adapting the defocus to the changing 
turbulence conditions, makes the final Strehl value increase to 0.45. The right panel of Fig. 4 
shows that the defocus increases when the seeing, and hence the phase distortions, increase. 
However, as the deformable mirror corrects the wavefront, the fluctuation sizes decrease and 
the defocus distance can again be reduced. This possibility of decreasing the defocus distance 
as the AO correction takes effect has been discussed F. Roddier [15]. 

In Fig. 5, the seeing varies randomly between 0.1" and 1.8". When the defocus is kept 
constant at its initial low value, the sensor measurements become dominated by diffraction 
as the seeing increases and erratic voltage commands distort the mirror surface. As the 
seeing decreases again, the distorted mirror surface still cause large diffraction patterns and 
the turbulence induced phase-distortions can not be measured: the AO loop has become 
unstable and the Strehl continuously decreases. When the defocus is adjusted, the average 
Strehl converges to the Strehl value (0.5) that roughly corresponds to the mean seeing (0.8"): 
indeed Fig. 2 shows that the Strehl equals 0.5 when the seeing equals 0.8" during the whole 
AO correction. 

Both simulations represent extreme cases, since the initial defocus lengths are obtained 
under excellent seeing conditions and are thus exceptionally small. The black curves of Fig. 6 
show that a defocus adjustment loop is not necessary for bright reference objects if one 
choses a constant and large defocus distance: the AO correction is then approximately the 
same for a constant and an adjusted defocus distance. Simulations and experimental results 
by O. Lai have likewise shown that in a low-order curvature system such as PUEO the AO 
loop performance is not very sensitive to changing defocus distances [16]. However, for faint 
reference stars it is not satisfactory anymore to keep the defocus distance at a constant large 
value (see Fig. 6): aliasing effects require the defocus distance to stay at its inferior limit, 
and an adjustment to the current seeing is required. This suggests that, as the number of 
correction elements of AO systems rise and each correction element has to work with less 
photons, defocus adjustment loops might become increasingly necessary. 

5. Conclusion 

We have tested two different methods to automatically adjust the defocus length of a cur- 
vature sensor to the continuously changing turbulence conditions. 

In the first method the length is adjusted prior to the AO loop through the acquisition 
of a long exposure image and is then kept constant. However in extended observations of 
astronomical objects or in any observation of fast moving targets, such as satellites, the 
turbulence intensity changes over the course of the AO sequence. For such observations, we 
have suggested a method to adjust the defocus length during the AO loop by use of the 
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voltage values sent to the deformable mirror. 

While the implications for AO in astronomy are apparent, it is also worthwhile in conclu- 
sion to note the potential interest of this work for ophthalmology. Detailed information on 
the retina requires observations from outside the eye. Without optical correction the quality 
of such observations is markedly degraded by aberrations due to the various dioptric surfaces, 
such as the cornea, and by the aqueous and the vitreous humour. The correction needs to 
be dynamic, because of the eye's motions and because the aberrations fluctuate. Hofer et 
al. [17] suggest a required AO correction speed of 10-40 Hz. 

Current methods to sense ocular aberrations, such as laser-ray tracing and spatially re- 
solved refractometry, are of limited applicability because they do not permit rapid acqui- 
sitions [18]. Dynamic aberration-correction is, therefore, only performed by use of Shack- 
Hartmann (S-H) sensors. Liang et al. [19] have shown that S-H sensors permit fast, precise 
and repeatable measurements of ocular aberrations. On the other hand, Moreno-Bruiso et 
al. [20] and Glanc et al. [21] found eyes where the S-H sub-aperture images are too distorted 
to be properly analyzed, and one needs to note that such eyes tend to be the most interesting 
from a clinical point of view. Porter et al. [22] characterized the distribution of the ocular 
aberration in a population of 109 subjects and they concluded, in line with previous studies, 
that the ocular aberrations differed greatly from person to person - even though their study 
did not include subjects with pathologies such as cataracts or keratoconus. This suggests 
that retinal imaging could benefit substantially from the use of AO systems in which the 
sensitivity can be easily adjusted to the eye's aberration and their temporal fluctuations. 

Curvature sensors have, in exploratory investigations, been successfully utilized to measure 
ocular aberrations [23], [24]. In these experiments, however, the sensors have not been used 
in combination with deformable mirrors, but were instead derived in terms of a complex 
wavefront reconstruction algorithm. Since curvature-based AO systems in astronomy have 
proven simple and efficient when used together with a bi-morphic or a membrane mirror [3], 
[5], [9], [10], [11], we suggest to implement a curvature based AO system for retinal imaging. 
Compared to a S-H based system, the chief advantage of such a system will be the capacity 
to rapidly adjust the defocus distance when different patients, with their particular ocular 
aberrations are examined. 
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Fig. 2. Strehl ratio and defocus distance versus seeing for a curvature sensor 
under various seeing conditions (See left column of Table 1 for the parameter 
values). Triangles: Sequences with different defocus values are simulated and 
the defocus that yields the highest Strehl is kept. Dashed line (resp. dotted 
and dashed- dotted): the minimum (resp. the average and maximum) defocus 
distance obtained with the former procedure is kept for all turbulence inten- 
sities. Solid line: the defocus value is determined from a long-exposure image, 
acquired prior to the AO loop. 
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Fig. 3. Seeing at 0.5/iin, recorded by the DIMM (dots) and MASS (squares) 
on the Mauna Kea observatory. The monitors are installed next to the 
Canada-France-Hawaii telescope which accommodates the curvature AO sys- 
tem PUEO. DIMM and MASS acquire a measurement about every two min- 
utes. MASS is insensitive to the lowest turbulent layers and generally measures 
smaller seeing values than DIMM. Credit: http://mkwc.ifa.hawaii.edu/ 
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Fig. 4. A 70- seconds AO sequence, during which the seeing, 6q, is incremented 
by 0.05" every 1000 th iteration. Left panel: Full line: Average Strehl when the 
defocus is adjusted during the sequence, dotted lines: the defocus is determined 
prior to the AO loop when 8q = 0.2". Circles: 8q/2. Right panel: Defocus 
values. Circles: 9 . (Parameter values: see right column of Table 1) 
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Fig. 5. A 70-seconds AO sequence, during which the seeing, 6 , is modified 
every 1000 th iteration. Left panel: Full line: Average Strehl when the defocus 
is adjusted during the sequence, dotted lines: the defocus is determined prior 
to the AO loop when 8q = 0.15". Circles: 9q/2. Right panel: Defocus values. 
Circles: 9 . (Parameter values: see right column of Table 1) 
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Fig. 6. A 70-seconds AO sequence, during which the seeing, #0, is modified 
every 1000 th iteration between 0.1" and 1.8". Full lines: Average Strehl when 
the defocus is adjusted during the sequence, dotted lines: the defocus is deter- 
mined prior to the AO loop when #0 = 1-6". The magnitude of the reference 
star varies between 5 (black), 9 (magenta) and 11 (blue). Circles: /2 (Fur- 
ther parameter values: see right column of Table 1). See the online edition of 
JOSA A for the color version of this figure. 
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Parameters 



Sec. 3 



Sec. 4 (if different from Sec. 3) 



Telescope diameter, D 
Diameter of central obstruction 
Size of influence functions, M x M 
Number of sub-apertures, N 

Size of phase screen 

Seeing at 0.5 /zm, # (0.5 /im) 

Layer speed, V 

Height of turbulent layer, H 

Outer scale of turbulence, L 

Target magnitude, m 
Sky magnitude, ms 
Photon flux at magnitude 
Atmospheric + instrumental throughput 
Imaging wavelength, Xj 
Sensing wavelength, \s 

AO loop frequency, / 

Read-out latency 

Read-out noise and dark current 

First ... iterations disregarded 

Minimal curvature radius 

of the deformable mirror 

Loop gain, g 

Extra-focal distance in an //60 beam, I 



7.9 m ; 220 pixels 
1.0 m ; 28 pixels 
237 x 237 pixels 
80 

2048 x 256 pixels 4096 x 256 pixels 

[0.2", 1.4"] [0.1", 1.8"] 

2m/s 
100 m 

+oo 

5 

20 

10 13 m -2 s -1 
0.25 
1.65 /im 
0.7 /mi 

1000 Hz 200 Hz 

0s 



50 

13 m 

[0, 1] increment: 0.1 

[0.1, 1.8] increment: 0.1 [0.1,3.0] increment: 0.1 



Table 1. Parameter values of the numerical simulations, if not otherwise spec- 
ified in the text. 
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